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Abstract 

The electron current density in nanoscale junctions is typically several 
orders of magnitude larger than the corresponding one in bulk electrodes. 
Consequently, the electron-electron scattering rate increases substantially 
in the junction. This leads to local electron heating of the underlying 
Fermi sea in analogy to the local ionic heating that is due to the increased 
electron-phonon scattering rates. We predict the bias dependence of local 
electron heating in quasi-ballistic nanoscale conductors, its effect on ionic 
heating, and discuss possible experimental tests of our results. 

In the process of electrical charge transport, the dissipation of energy via 
heat production plays a significant role. This effect is of particular importance 
in nanoscale systems, like e.g., atomic or molecular structures between bulk 
electrodes, ^ since it determines their structural stability under current flow. It 
is now understood that the large current densities in nanojunctions, compared 
to their bulk counterparts, may lead to substantial heating of the nanostructure 
ions.— 1^1^ This effect is directly related to the consequent increase of the electron- 
phonon scattering rates in the junction. In quasi-ballistic systems, i.e. when 
the mean free path is much longer than the dimensions of the nanostructure, by 
assuming a bulk lattice heat conduction mechanism, the local ionic temperature 
is predicted to be Tion oc \/V, where V is the applied bias.— i^i^ For the same rea- 
sons, and due to the viscous nature of the electron liquid,— we here suggest that 
the local increase of the electron-electron scattering rate in the junction gives 
rise to local heating of the underlying Fermi sea, whether the system has one or 
many conducting channels. This local electronic temperature would also affect 
the electron-phonon scattering rates, and consequently the bias dependence of 
the ionic temperature. 

In this Letter, we first estimate the bias dependence of the local electron 
temperature in quasi-ballistic systems, assuming no ionic heating is produced. 
We will then determine the effect of the electron heating on the bias dependence 
of the ionic temperature. We finally discuss possible experiments to test our 
predictions. 
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Let us start from a simple argument on the expected bias dependence of the 
electron heating. Assume first no electron-phonon scattering is present in the 
system. The power generated in the nanostructure due to exchange of energy 
between the current-carrying electrons and the underlying Fermi sea has to be a 
small fraction of the total power of the circuit V"^ /R {R is the resistance). Let us 
define this fraction as P = aV^/R, with a a positive constant to be determined 
from a microscopic theory. At steady state this power has to balance the thermal 
current Ith (heat per unit time) carried away into the bulk electrodes by the 
electrons. Let us first assume that the electron thermal conductivity k follows 
a bulk law, i.e. it is related to the specific heat per unit volume cy via the 
relation k — ypXeCy/d where vp is the Fermi velocity, Ae is the electron mean 
free path, and d the dimensionality of the system.— We will later give a general 
argument to justify the form of this law in quasi-ballistic systems away from 
electronic resonances. We know that at small temperatures the specific heat of 
an electron liquid is proportional to the electronic temperature Te so that 

k = iT, (1) 

where 7 = kpkgXe/9h in 3D, and 7 = TrkpkgXe/Gh in 2D, kp is the Fermi 
momentum, and ks the Boltzmann constant . The thermal current, Ith oc fcTg, 
is then given by Ith = I'T"^ where 7' has to be determined from a microscopic 
theory. At steady state the condition P = Ith implies a linear increase of the 
electronic temperature with bias, 

Te = 7e-eV, (2) 

if the coefficient je-e does not vary appreciably with bias. In macroscopic 
electrical contacts this result is known as the 'V — Q relation" .— Here, we will 
re-derive it from a microscopic theory and determine the quantity 7e-e in terms 
of the properties of the junction and of the electron liquid. This theory needs to 
take into account the influence of electron-electron interactions on both the heat 
production and dissipation. Two of the authors (RD'A and MDV) have recently 
shown that the dynamics of the electron liquid in nanojunctions can be described 
using a hydrodynamic approach,— (see Ref. (9) for a different hydrodynamical 
approach to the electron liquid flow) so that the equations of motion for the 
electron density n(r,t) and velocity field v(r,t) can be written in the Navier- 
Stokes form^iia 

Dtn{r, t) = - n(r, t) V • v{r, t) 
mn{r,t)DtV,{r,t) ^ - V ,P{r,t) +V j(T,j{r,t) (3) 
-n{r,t)'sJ,Ve^t{r,t). 

where P is the pressure, Dt — dt + v ■ V \s the convective derivative operator, 
and Vext is an external potential (like the electron-ion potential). The stress 
tensor (7^ ^ is — i^iiiii^ 

o-ij = ri (diVj{r,t) + djV^{r,t) ~ • v \ + C%V • u, (4) 
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where rj and ^ are positive parameters which correspond to the shear and bulk 
viscosity of the hquid, respectively. They have been determined for the electron 
liquid using linear response theory.— Since in the DC limit (the regime of interest 
here) C <C in what follows we will consider the shear viscosity only. This term 
is the one responsible for the heating of the Fermi sea: electrons crossing the 
junction experience an internal "friction" with the electrons of the underlying 
Fermi gas, thus creating electron-hole pairs and generating heat. In keeping 
with this hydrodynamic picture we supplement Eqs.[n]for the particle densities 
and velocity field with an equation which describes the energy dissipation and 
diffusion. Since the quantity Vidi j- gives the force acting on the electron liquid 
caused by viscosity, the equation for the heat transfer read&^ 

Te{r,t)Dtsir,t) =a^jir,t)djVi{r,t) 

+ V- [fc(r,i)Vr,(r,i)] ^' 

where k is again the thermal conductivity and s is the local entropy per unit 
volume. 

For metallic quantum point contacts (QPC) of interest here the electron 
liquid can be assumed incompressible. ^'^ At steady state eq. [5] therefore becomes 

a,,,(r)aj«,(r) + V • [fc(r)Vr,(r)] = cv{Te)v{r) ■ VT,{r). (6) 

where we used T{r)\I s{r) ~ cp\IT{r), Cp being the specific heat per unit volume 
at constant pressure.— For an electron liquid at low temperatures cy ~ cp.— 
Equations [3] and El with all quantities time- independent, together with the def- 
inition of the stress tensor [4l constitute a set of equations that describes the 
charge and heat flow at steady state. In this effective theory quantum mechan- 
ics enters explicitly in the electron liquid constants rj, k, and cy.^^ 

Before proceeding with the calculation of the temperature from eq. [S] we 
need to know the form of the thermal conductivity k. We give here a simple ar- 
gument to justify the use of the bulk form, eq.[T]in quasi-ballistic systems at low 
temperatures (and bias) and in the absence of electronic resonances. Following 
the Landauer approach,- - let us assume the electrons tunnel across the junction 
from the right electrode which is in local equilibrium with temperature Tr and 
chemical potential /i_R, to the left electrode in local equilibrium with tempera- 
ture Ti and chemical potential /i^,. If the transmission coefficient for electrons 
to move elastically from right to left is T{E), the energy current transported 
by the charge is of the form 

Ith ^ JdE ET{E)[fR{E, Tr) - fUE, Tl)] (7) 

where fR[L) is the local Fermi-Dirac distribution function of the right (left) 
electrode at its own local temperature and chemical potential. Let us set = 
and consider a small right electrode temperature Tr = T^. For small biases, 
in the absence of electronic resonances, the transmission coefficient T{E) can 
be assumed independent of energy, T{E) = T . From eq. [7| we therefore see 
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Figure 1: Schematic of the local electronic temperature profile of a nanojunction. 

that Ith T^, i.e. the thermal conductivity must be of the form [TJ In the 
following we will assume T ~ 1, typical of metallic QPCs, and therefore assume 
the form [1] for k with the same coefficients. 

In order to obtain the temperature profile of the electron gas inside the 
nanojunction and evaluate the relation between and V , we need to solve eq. 
[Altogether with Navier-Stokes equations[3]with the boundary conditions imposed 
by the geometry of the junction. For an arbitrary system, this is obviously an 
impossible task. Since we are only interested in the analytical dependence of 
the temperature on bias, and the estimate of the coefficients in this relation, we 
proceed as follows. Let us assume a given profile for the velocity of the fluid 
and solve eq. [6] with the boundary condition that the temperature reaches the 
constant bulk value far away from the constriction. We choose a reference frame 
in which the direction of the current flow is along the x axis and the constriction, 
whose length is L, is centered at a; = (see Fig. [T|). Let us also assume that the 
velocity varies negligibly in the (y, z) directions compared to the much larger 
variation along the x direction. 

Since the system has a variable cross-section and the total current / has to 
be constant throughout the device, the fluid velocity and the cross-section of 
the system are approximately related via 



The cross-section profile A{x) contains the information about the geometry of 
the device. In the following we assume the simplest profile possible which allows 



/ = envx{x)A{x). 



(8) 
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an analytical solution: the adiabatic model (see Fig. [T|) 



A. + r- 



- ^ ) (9) 



with Ac and r positive parameters describing the minimal cross-section of the 
junction and the "rate" at which it opens, respectively, and 6{x) is the Heaviside 
step function. With this choice of the profile A{x) we get the velocity profile 
we need to supply in eq. [S] to obtain the temperature profile 
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The total electron velocity u{x) is the sum of the Fermi velocity and the 
fluid velocity, u = vp + Vx- In the bulk where the heating can be considered 
negligible, u is related to the applied bias voltage via mw^ = 2{Ep + eV), and 
since in general, <C vp we get Vx — eV/mvp for the fluid velocity in the 
bulk.— This is the boundary condition that eq. [10] has to satisfy for |x| > L/2, 
i.e. Ac + rL'^/A = mvpl /ne^V . This condition implies rh"^ /AAc ^ 1 which 
simply reflects the physical fact that, due to screening, the electron velocity 
approaches the bulk value very fast away from the junction. 

In the following we assume that deep inside the electrodes the electron gas is 
in equilibrium with a zero-temperature thermal bath, i.e. lim|2.|^oo Tc{x) = 0.— 
The solution to eq. [S] can be calculated analytically and expressed in terms of 
a combination of rational, trigonometric and hypergeometric functions. From 
this solution we can estimate that the contribution of the the advection term 
cvv-VT is controlled by the dimensionless quantity F = {Id/vp^eneAc)^ Ac/r. 
Due to the quasi-ballistic assumption, we find F <C 1 so that the advection term 
can be neglected. 

In the absence of the advection term the solution to eq. [S] inside the con- 
striction (|a;| < L/2) is given by 
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Let us then estimate the maximum temperature which, in the absence of the 
advection term, occurs at x = 0, 




— 1 77 / / r L , 
3d j\ y Ac2 
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(1 +rL2/4Ac)2 

(12) 

In the linear transport regime and for T ~ 1, / = GqV with Gq = 2e^ jh. As 
we have anticipated in eq. [H the relation between the temperature and the bias 
then simplifies to Tm = 7e_e^, where 7e-e is a constant that can be read from 
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eq. [T21 It is interesting to note that a similar expression has been found in the 
case of superconducting junctions, where, however, the maximum temperature 
inside the junction is found to be independent of both the junction geometry 
and the electron properties.— 

By using the expression of 77 as a function of the particle density given in 
Ref. (5) we can estimate the value of the constant 7e-e for various systems. We 
consider here both a 3D gold point contact and a 2D electron gas (2DEG). For 
a 3D gold {vs = 3) QPC with effective cross-section Ac = 7.0 A^, by assuming 
a typical inelastic mean free path in quasi-ballistic systems of Ag ~ 100 nm, we 
get 7e-e(QPC) ~ 65 K/V. For a 2DEG, assuming — 10, Ae ~ 10 /^m, and 
Ac = 20 nm we get jc-ei^DEG) ~ 1.2 x 10^ K/V. 

Let us now discuss the effect of local electron heating on the ionic heating. In 
the 2DEG we expect a negligible heating of the ions. In atomic-scale junctions, 
instead, the ions may heat up substantially due to electron-phonon scattering. 
As discussed at the beginning of the paper, the ionic temperature in quasi- 
ballistic systems is predicted to be Tion = "Je-pVV, if no electron heating is 
taken into account. The constant 7e-p has been estimated in— for atomic and 
molecular systems. If we now allow for the electron heating, then part of the 
ionic heating is lost in favor of the local electron temperature. Since the source 
of energy is the same for both processes (the kinetic energy of the electron), the 
ionic temperature must be lower in the presence of the electron heating than 
without it. To first order we can assume the electron-electron and electron- 
phonon processes independent and occurring with equal probability. Let us also 
suppose that, for a given bias, the local electron temperature in the absence 
of electron-phonon scattering, Tg, is smaller than the ionic temperature in the 
absence of electron heating, Tion- The energy lost by an electron due to heating 
is "seen" by the phonons as a "sink" of energy at the electron temperature Tg. 
This energy however is "lost" by the phonons at their own rate. The balance 
between the power absorbed by the phonons in the junction and the power 
dissipated away at steady state therefore gives the new ionic temperature 

T^^[jt^V'-jlcV'f'. (13) 

If Tion — Tg the processes of energy release from current-carrying electrons to 
the phonon gas or to the Fermi sea cannot be considered independent and a 
microscopic theory that considers them on equal footing is needed. 

We point out that for metallic QPCs the theoretical values for je-p are of 
the order of hundreds of K/V^/^ while our estimate for je-e is of the order of 
tens of K/V.^'- We thus expect that, for these systems, the electron heating 
contribution to the ionic temperature is small for a wide range of biases. In the 
case of a 2DEG, on the other hand, the ionic heating is either not present or it is 
very small and it may be possible to measure the effects of the electron temper- 
ature directly. In the case of organic molecules between metallic electrodes it is 
known that je-p decreases exponentially with the length of the molecules.'^ For 
such systems, we thus expect the two effects to be of similar importance and one 
should therefore be able to observe deviations from the ^/V dependence of the 
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ionic temperature due to electron heating as predicted in eq. [T31 However, the 
exact value of the parameter 7e-e requires generalization of the present theory 
to the case of non-ideal conductance, i.e., for a transmission coefficient T less 
than 1. This can, in principle, be done but no general analytical solution can be 
derived. The reason is that one needs to know both the thermal conductivity 
dependence on the transmission coefficient, and the junction profile that leads to 
such coefficient. For T <C 1 this may even lead to a non-linear dependence of Tf. 
as a function of V . We can however comment that if one works at fixed current 
while varying T, for example, by acting with a gate voltage or by stretching 
the chemical bonds of the molecule, an increase of electron temperature due 
to reduced thermal conductivity is expected. We finally suggest that the local 
increase of the electron temperature in a nanojunction may be extracted by 
measuring the Johnson-Nyquist noise.— Any deviation of this noise compared 
to the expected one in the absence of electron heating may give a direct evidence 
of this effect. 

In conclusion, we have discussed the phenomenon of local electron heating 
in nanoscale junctions. We have predicted its bias dependence and estimated, 
using a hydrodynamic approach, its magnitude in quasi-ballistic systems. We 
have also discussed its role on the local ionic heating and suggested systems 
where this effect is most likely to be observed. 
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